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Abstract: Recently the gaussian expansion method has been apphed to investigate the 
dynamical generation of 4d space-time in the IIB matrix model, which is a conjectured 
nonperturbative definition of type IIB superstring theory in 10 dimensions. Evidence for 
such a phenomenon, which is associated with the spontaneous breaking of the SO(IO) 
symmetry down to S0(4), has been obtained up to the 7-th order calculations. Here we 
apply the same method to a simplified model, which is expected to exhibit an analogous 
spontaneous symmetry breaking via the same mechanism as conjectured for the IIB matrix 
model. The results up to the 9-th order demonstrate a clear convergence, which allows us 
to unambiguously identify the actual symmetry breaking pattern by comparing the free 
energy of possible vacua and to calculate the extent of "space-time" in each direction. 
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1. Introduction 

It has been long considered that matrix models may be useful as a nonperturbative formu- 
lation of string theory, and hence play an important role similar to the lattice formulation 
in quantum field theory. Indeed matrix models have been quite successful in formulating 
non-critical string theory, and after the development of the notions such as "string duality" 
and "D-branes" , the idea has been extended also to critical strings. The IIB matrix model 
[1] is one of such proposals, which is conjectured to be a nonperturbative definition of type 
IIB superstring theory in 10 dimensions. It is a supersymmetric matrix model, which can 
be formally obtained by the zero-volume limit of lOd SU(A'^) super Yang-Mills theory. 

In this model the space-time is represented by the eigenvalue distribution of ten bosonic 
matrices [2]. If the distribution collapses dynamically to a four-dimensional hypersurface, 
which in particular requires the SO (10) symmetry of the model to be spontaneously broken, 
we may naturally understand the dimensionality of our space-time as a result of the nonper- 
turbative dynamics of superstring theory. In ref. [3] the first evidence for the above scenario 
has been obtained by calculating the free energy of space-time with various dimensionality 
using the gaussian expansion method up to the 3rd order. Higher-order calculations [4, 5] as 
well as the tests of the method itself in simpler models [6, 7] have strengthened the conclu- 
sion considerably. Refs. [8] provide another evidence for the emergence of four-dimensional 
space-time based on perturbative calculations around fuzzy-sphere like solutions. 

While these results are certainly encouraging, it is desirable to understand the mecha- 
nism for the spontaneous symmetry breaking (SSB) of the rotational symmetry. In refs. [9] 
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it has been pointed out that the phase of the fermion determinant favors lower dimensional 
configurations since the phase becomes stationary around such configurations ^. Indeed 
Monte Carlo simulations show that SSB does not occur in various models without such a 
phase factor [12, 13, 14]. Unfortunately including the effects of the phase in Monte Carlo 
simulation is technically difficult due to the so-called sign problem, but a new method 
[15], which is tested in Random Matrix Theory [16], was able to produce some preliminary 
results, which look promising. In ref. [17] a simple matrix model which realizes the above 
mechanism has been proposed. The model contains A'f flavors of Weyl fermion in the fun- 
damental representation of SU(A^), which yield a complex fermion determinant, and the 
large N limit is taken with the ratio r = Nf/N being fixed. The model can be solved 
exactly at infinitesimally small r, and the SO (4) symmetry is shown to be broken down to 



In this paper we study this model at finite r by the gaussian expansion method. Since 
the model is much simpler than the IIB matrix model, we can perform calculations up to 
the 9-th order with reasonable efforts. The results demonstrate a clear convergence for 
r < 2, which allows us to unambiguously identify the symmetry breaking pattern and to 
calculate the extent of "space-time" in each direction. 

In fact it turns out that the S0(4) symmetry is broken down to S0(2) at finite r. 
However, at small r we reproduce the free energy as well as the extent of "space-time" 
in each direction obtained in ref. [17], which implies that the S0(3) symmetry is realized 
asymptotically as r approaches zero. In the large r region, on the other hand, the extent 
of "space-time" in two directions, in which the SO (2) symmetry is realized, becomes much 
larger than the remaining two directions. This behavior can be understood from the 
viewpoint of refs. [9] since the phase of the fermion determinant becomes stationary for 
two-dimensional configurations, and increasing r tends to amplify the effect of the phase. 
Thus our results nicely demonstrate the proposed mechanism for the dynamical generation 
of space-time in the IIB matrix model. 

The rest of this paper is organized as follows. In Section 2 we define the model and 
review the known results. In Section 3 we explain how to apply the gaussian expansion 
method to the model. In Section 4 we present our results. Section 5 is devoted to a 
summary and discussions. The details of our calculations are given in the Appendix. 

2. The model 

The model we study in this paper is defined by the partition function [17] 



S0(3). 




(2.1) 



(2.3) 



(2.2) 



^See refs. [2, 10, 11] for discussions on other possible mechanisms. 
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where Af^ = 1, • • • ,4) arc N x N tracelcss ^ hermitian matrices and ipa, ipa (ct = 1,2; 
/ = !,•••, A'^f) are A^-dimensional row and column vectors, respectively, making the system 
SU(iV) invariant. The integration measure for is given by 

a=l /i=l 

where A^^^ is the coefficient in the expansion Af, = Y.a=i^ K ^i*^ respect to the SU(iV) 
generators (a = 1, • • • , (iV^ — 1)) normahzed as tr (T'^T'') = ^S"-^ . The integration 
measure for the fermions is given by 

Nf N 2 

dV'dvi=nnnd^« . (2.5) 

/=1 i=l a=l 



The system has an S0(4) symmetry, under which transforms as a vector, and ipa 
and transform as Weyl spinors. The 2x2 matrices T ^ are the gamma matrices after 
the Weyl projection, and an explicit form is given for instance by 

(2.6) 






The fermionic part of the model can be thought of as the zero- volume limit of the system of 
Weyl fermions in four dimensions interacting with a background gauge field via fundamental 
coupling. 

We take the large limit keeping the ratio r = Nf/N fixed (Veneziano limit). In 
order to discuss the SSB of the S0(4) symmetry in that limit, we consider the "moment 
of inertia tensor" [2, 12] 

r^. = ^tr(^^^,) , (2.7) 

which is a 4 X 4 real symmetric tensor, and denote its eigenvalues as {A^ ; i = 1, • • • , 4} with 
the specified order 

Ai > A2 > A3 > A4 . (2.8) 

If the vacuum expectation values (VEVs) (Aj) (i = 1, • • • , 4) do not agree in the large N 
limit, we may conclude that the SSB occurs. Thus (Aj) plays the role of an order parameter. 
In the present model the sum of the VEVs is given by 

i=l //=1 ^ / \ / 

for arbitrary N and r due to a "virial theorem" [17]. 



^The tracelessness condition was not imposed in the original paper [17]. While this condition does not 
affect the large N limit of the model, it simplifies our calculation drastically; see footnote 4. 
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At infinitesimally small r the VEVs can be obtained in the large N limit as [17] 

(Ai) = (A2) = (A3) = l + r + o(r) , 

(A4) = l-r + o(r) , (2.10) 

which means that the SO (4) symmetry is spontaneously broken down to SO (3). The SSB 
is associated with the formation of a condensate (ipitpa), which is invariant under S0(3), 
but not under the full S0(4) transformation. 

An important feature of the model that is relevant to the SSB is that the fermion 
determinant detP, where P is a 2N x 2N matrix given by P = T^^d^, is complex in 
general. If one replaces the fermion determinant by its absolute value, the same analysis 
at infinitesimal r leads to an S0(4) symmetric result [17]. Thus the SSB of the original 
model occurs precisely due to the phase of the fermion determinant. 

At large r the effect of the phase is amplified, and we may expect that the configurations 
for which the phase is stationary dominate the path integral. Analogously to the situation 
in the JIB matrix model [9], the phase becomes stationary for 2-dimensional configurations 
in the present model. Therefore we anticipate the emergence of 2-dimensional "space-time" 
(Ai, A2 3> A3, A4) as r increases ^. 

3. The gaussian expansion method 

We are going to obtain results for the model (2.1) at finite r using the gaussian expansion 
method. The method has a long history, and the original idea appeared already around 
1980 in the context of solving quantum mechanical systems [18, 19], where the expansion 
was shown to be convergent in some concrete examples [20]. The method proved useful 
also in field theories [21] in various contexts. Applications to supcrstring/M theories using 
their matrix model formulations have been advocated by Kabat and Lifschytz [22], and the 
subsequent scries of papers [23] revealed interesting blackhole thermodynamics. Applica- 
tions to simplified versions of the IIB matrix model were initiated in refs. [24]. An earlier 
application to random matrix models can be found in rcf. [25]. 

Similarly to the case of the IIB matrix model [3, 4, 5], let us introduce the gaussian 
action ^ 

4 N( 2 

So = ^NY,t,triA,f + NY, E ^-p^i^. (3.1) 

M=l /=1 a,l3=l 

which breaks the S0(4) symmetry. The 2x2 complex matrix A can be expanded in terms 
of gamma matrices as 

4 
/i=i 

^Note that the phase of the fermion determinant is invariant under the scale transformation Ay, 1— > aAij,. 
Therefore it is only the ratio of the eigenvalues that matters for the stationarity. 

* If we did not impose the tracelessness condition on A^, we would need to consider a linear term such 
as Siin = ^"12 fx hfitr Afx in the gaussian action (3.1). The gaussian expansion method can be extended to 
such a case [7], but the calculation will be more involved. 
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using 4 complex parameters U/j,. Then we consider the action 

SGEM{t,u;X) = J [[So + A(5b - So)] + 5f] , (3.3) 

which reduces to the original action for A = 1. The gaussian expansion amounts to cal- 
culating various quantities as an expansion with respect to A up to some finite order and 
setting A = 1 eventually. As we will discuss shortly, the free parameters t/^ and in the 
gaussian action play a crucial role in the method. 

In fact the gaussian expansion can be viewed as a loop expansion with the "classical 
action" (Sq + Sf) and the "one- loop counterterms" (S'b — Sq)- This becomes clear upon 
rescaling A^j^ and as i— > A^^ tpa '—^ V^i/Ja, so that the partition function takes the 
form 

Z = J dAdipdipe-^^'^'+^'-'\ (3.4) 

S,i{t, u) = So + VXSi , 5e.t.(t, ^^) = A (Sb - 5o) . (3.5) 

In actual calculations the "one-loop counter terms" can be incorporated easily by exploiting 
the relation 

Sci{t,u) + Sc.tXt,u) = Sd{t + X{1 - t),u - Xu) . (3.6) 

As an example, let us consider evaluating the free energy F = -j^lnZbythe gaussian 
expansion method. We first calculate the free energy J^(t, u) for the "classical action" 
Sc\it,u) defined by 

exp[-N^Tit,u)] = J dAdV'dVie-^-i^*'") . (3.7) 

This can be done by ordinary Feynman diagrammatic calculations, where the use of 
Schwinger-Dyson (SD) equations reduces the number of diagrams considerably [4]. Suppose 
we obtain the result up to the K-th order as 

K 

TK{t,u) = Y,^k{t,u)XK (3.8) 

fc=0 

We shift the arguments, and obtain the new coefficients Fk{t,u) in the expansion 

K 

J^K{t + A(l -t),u- An) = Fk{t, u) X'' + 0{X^+^) . (3.9) 

A;=0 

Then the free energy for the original model can be evaluated as 

K 

FK(t,u) = Y,Fkit,u) . (3.10) 

fc=0 

The result of such calculations depends on the free parameters in the gaussian action. 
However, in various models [4, 6, 7] a "plateau" region, in which the result becomes almost 
constant, was found to develop in the parameter space as one goes to higher orders of the 
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expansion. Moreover it turned out that the height of the plateau agrees very accurately 
with the correct value obtained by some other method. Therefore it is reasonable to expect 
that the method works in general if one can identify a plateau in the parameter space. In 
old literature the free parameters were determined in such a way that the result becomes 
most insensitive to the change of the parameters [19], but it is really the formation of a 
plateau that ensures the validity of the method as has been first recognized in ref. [4]. 

Identification of a plateau becomes a non-trivial issue when there are many parameters 
in the gaussian action. The histogram prescription [6, 7], which works nicely when there 
axe only one or two real parameters, does not seem to work when the number exceeds three. 
(Note that we have 4 real and 4 complex parameters in the present case.) We have also 
attempted a Monte Carlo simulation in the parameter space to search for a plateau but 
with little success. We therefore use the prescription adopted for the IIB matrix model. 
First we solve the "self-consistency equations" 

|;F;,(t,.)=0, 

^F,(t,„)=0. (3.11) 

Typically we obtain many solutions as we go to higher orders. If we observe that solutions 
concentrate in some region of the parameter space, we consider it as an indication of the 
plateau formation. 

Although the number of parameters is much less than that (10 real and 120 complex 
parameters) in the IIB matrix model, it is still difficult to obtain all the solutions of the 
self-consistency equations at high orders. As is done in the case of the IIB matrix model [3, 
4, 5], we search for solutions assuming that some subgroup of the full rotational symmetry 
is preserved. Here we consider the following Ansatze. For each case the independent 
parameters will be 2 real and 1 complex numbers. 

SO(3) Ansatz : We assume S0(3) symmetry in the X2, 3:3, directions. Then the 
parameters are restricted to 

^2 = ^3 = i4(= i), U2 =U3 = U4 = . (3-12) 

SO(2) Ansatz : We assume S0(2) symmetry in the X3, X4 directions. Furthermore we 
impose discrete symmetry under xi X2, X2 xi, X4 —X4. Then the parameters are 
restricted to 

tl=t2, ts=t4(=i), Ui=U2, Us = U4 = 0. (3.13) 



4. Results 

For each solution of the self-consistency equations obtained within the symmetry Ansatz, 
we calculate the free energy. (See appendix for the details.) The free energy we plot in 
what follows is actually defined by 

/= lim {F-2(l-r)lniV} , (4.1) 

W— ►oo 
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where the subtraction is necessary to make the quantity finite. The exact result at in- 
finitesimal r is given by [17] 

/ = -21n2 + (l-ln2)r + o(r) . (4.2) 




0.5 1 1.5 2 0.5 1 1.5 2 




0.5 1 1.5 2 0.5 1 1.5 2 



r r 

Figure 1: The free energy obtained for the SO (3) Ansatz at orders 3(left top), 5(right top), 7(left 
bottom) and 9(right bottom) is plotted as a function of r. 

In Fig. 1 the free energy calculated for the SO (3) Ansatz at orders 3,5,7,9 is plotted 
against r. We find that at orders 5 and 7 there are two solutions which almost coincide 
with each other throughout the whole region ^ of < r < 2. At the 9-th order there 
are actually three solutions lying on top of each other, which are represented by the solid 
lines to be distinguished from the other solutions. We consider this as an indication of the 
plateau formation. Similar behavior is observed for the S0(2) Ansatz as one can see from 
Fig. 2. Again the three solid lines in the right bottom plot represent the solutions that we 
consider to be concentrating. 

We pick up the three solutions which concentrate at the 9-th order for the two Ansatze 
and plot them in Fig. 3 for comparison. Throughout the whole range of r considered, 
the free energy for the S0(2) Ansatz is smaller than that for the S0(3) Ansatz. Thus we 

^At r > 2 the solutions that concentrate at the 9-th order start to separate, and we cannot obtain 
reliable results in that regime. Note in this regard that the number of fermion loops gives the power of r 
in the result of the gaussian expansion. Therefore it is reasonable that the convergence becomes slower as 
we go to larger r. 
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conclude that the true vacuum is described by the SO (2) Ansatz. On the other hand, the 
results for the two Ansatze asymptote to each other as r approaches zero. The meaning of 
this behavior will be clarified shortly. 



-1 
-1.5 
-2 

^ -2.5 
-3 
-3.5 
-4 

0.5 1 1.5 2 

r 

Figure 3: Comparison of the free energy obtained for the SO (2) and SO (3) Ansatze (solid lines 
and dashed lines, respectively). The solutions that concentrate at the 9-th order are extracted from 
Figs. 1 and 2. 
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Let us move on to the calculation of observables. Similarly to the free energy, we can 
calculate an observable as an expansion with respect to A using the action (3.3). In our 
model the observable of primary interest is the eigenvalues Aj of the "moment of inertia 
tensor" (2.7). We calculate them for the S0(3) and S0(2) Ansatze at the 9-th order as 
a function of the free parameters in the gaussian action, and plug in the three solutions 
that are seen to concentrate in the study of free energy. In fact the VEVs of A^ can be 
readily obtained by diagonalizing Cfj^u defined by eq. (A.l), which is calculated anyway in 
the calculation of free energy ^. Fig. 4 shows the results. Note that we are ultimately 
interested in the results for the SO (2) Ansatz since it gives the smaller free energy. 

For the S0(3) (S0(2)) Ansatz the lines that grow almost linearly actually represent 
three (two) eigenvalues, which are degenerate due to the assumed symmetry. For the SO (2) 
Ansatz it turns out that the third largest eigenvalue comes closer to the two degenerate 
largest ones as one approaches r = 0, thus realizing the S0(3) symmetry asymptotically. 
In fact the results obtained for the SO (2) Ansatz are indistinguishable from those obtained 
for the S0(3) Ansatz at small r. This is consistent with our observation in Fig. 3 that the 
free energy for the S0(2) Ansatz have the same asymptotic behavior for r — as that for 
the S0(3) Ansatz. Actually we find that both the free energy and the observable agree 
asymptotically with the exact results (4.2), (2.10) at infinitesimal r. 

At large r, on the other hand, the results for the SO (2) Ansatz show a clear tendency 
that the two degenerate eigenvalues become much larger than the other two, which implies 
the emergence of a two-dimensional "space-time" . This agrees with the argument based on 
the phase stationarity given at the end of Section 2. Thus the gaussian expansion method 
enables us to obtain results at finite r, which naturally interpolate the SO (3) symmetric 
exact result at infinitesimal r and the two-dimensional behavior expected at large r. 
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Figure 4: The four eigenvalues of the "moment of inertia tensor" (2.7) obtained for the S0(3) 
Ansatz (left) and the S0(2) Ansatz (right) at the 9-th order are shown as a function of r. Note 
that the largest eigenvalue has 3-fold (2-fold) degeneracy for the SO (3) Ansatz and the SO (2) 
Ansatz, respectively. The three types of line correspond to the three solutions of the self- consistency 
equations (3.11) that concentrate at the 9-th order. 



^Sincc c^,y takes the form (A. 17) and (A. 20) respectively for the S0(3) and SO(2) Ansatze, the diago- 
nalization is actually needed only for the SO(2) Ansatz. 
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5. Summciry and Discussions 

In this paper we have apphed the gaussian expansion method to a matrix model which is 
expected to exhibit SSB of rotational symmetry due to the phase of the fermion determi- 
nant. The free energy calculated by the gaussian expansion method depends on the free 
parameters in the gaussian action, and the formation of a plateau in the parameter space is 
crucial for the validity of the method. For each Ansatz considered for the possible breaking 
pattern of the SO (4) symmetry, we obtained a clear evidence for the plateau formation. By 
comparing the free energy obtained for each Ansatz, we concluded that the true vacuum is 
described by the S0(2) Ansatz. Our results for the extent of "space-time" in each direction 
arc consistent with the exact result for infinitesimal r and also with the behavior expected 
at large r. 

The mechanism for the SSB of rotational symmetry demonstrated in the present model 
is expected to be at work also in the IIB matrix model. However, we should also note the 
difference of the two models. In the present model the fermionic degrees of freedom match 
with the bosonic ones at r = 1, but there the "space-time" is not really two-dimensional, 
but it looks like a "rugby ball" with two directions more extended than the other two. 
Moreover, if we consider the ten-dimensional version of the present model, we expect that 
there will be eight directions more extended than the other two. In the IIB matrix model, 
on the other hand, the result of the gaussian expansion method shows that the ratio of the 
extent of space-time in four directions to that in the remaining six directions increases with 
the order up to the 7-th order [3, 4, 5]. This suggests that four directions are much more 
(possibly, infinitely more) extended than the remaining six directions. It is conceivable that 
supersymmetry plays an important role here. Let us recall that the effective theory for the 
eigenvalues of the bosonic matrices in the IIB matrix model is a weakly bound system like 
a branched polymer due to cancellation between the bosonic and fermionic contributions 
[2]. This makes the space-time easy to collapse. Note also that the convergence of the 
gaussian expansion is not so clear in the IIB matrix model as in the present model. 

While it is certainly worth while to proceed to the 8-th or 9-th orders in the IIB matrix 
model, we consider that Monte Carlo simulations along the line of ref. [15] are necessary to 
definitely confirm the emergence of a four-dimensional space time. That approach is also 
expected to provide an intuitive understanding of why "4" instead of 3 or 5. Note in this 
regard that the present simplified model shares the technical difficulty for implementing 
the phase of the fermion determinant in Monte Carlo simulation. Our new results obtained 
in this paper therefore provide a nice testing ground for the new method for simulating the 
IIB matrix model. 
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A. Details of the calculation 

As explained in Section 3 the main part of calculations in the gaussian expansion is actually 
nothing but the ordinary perturbative calculation of the free energy T{t,u) defined by 
(3.7). The Feynman rules are given in Fig. 5. The first and second lines represent the bare 
propagators {{Af^)ij{A^)ki)Q and (V'a*^|"')oi respectively, where the symbol ( • )o represents 
a VEV obtained with the gaussian action (3.1). The third line stands for the interaction 
vertex coming from Sf. Instead of evaluating all the diagrams contributing to J^{t,u), we 
use the SD equations to reduce the number of the diagrams to be computed following the 
idea put forward in ref. [4]. 

i 1 

IIIIIIIIIIv 
j k 



( Si 



N 



1 1 

t ., 



H V 



— §.. 5f„ (A 



1/2 

X N ( r ) 



ij kl fg 



Figure 5: The Feynman rules for the A-expansion of (3.7). 



Note that the one-point function {{A^)ij) is proportional to Sij {tr A^) because of the 
SU(A'^) symmetry and therefore it vanishes due to the tracelessness condition imposed on 
A^. The full propagators (connected two-point functions) can be written in the form 

{{An)ij{A^)ki) = — (^SiiSjk - j^SijSki^ Cfj_u , (A.l) 



- V + 0(A) , V = A-' + 0(A) . (A.3) 



where the coefficients c^^ and V are given at the leading order in A as 

1 

Corresponding to eq. (3.2), the full propagator V can be parametrized as 

4 

V = ^di,f^, (A.4) 
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where the "conjugate" gamma matrices are defined by 

ffc = rfe (fc = i,2,3), r4 = -r4. (A.5) 

In what follows, we restrict ourselves to the large-A^ limit with the ratio r = ^ fixed, so 
that we have to consider planar diagrams only, but the method itself is applicable to finite 
as well. 

By using the SD equations for the full propagators, we can reduce the calculation of 
the free energy to that of two-particle-irreducible (2PI) planar diagrams. Here, by "2PI 
diagrams" we mean those diagrams which cannot be separated into two disconnected parts 
by cutting two propagators. Let us consider the 2PI planar vacuum diagrams whose internal 
lines are all replaced by the full propagators. The sum of such diagrams is a function of 
c^i, and and shall be denoted as N'^G{c, d). For example, N'^G{c, d) up to the 5th order 
can be computed from the 5 diagrams in Fig. 6, and it is given explicitly as 

G(c, d) = - ^ A r ^ c^, IV (f^PF.p) 

+ ^^^^' c^^capC«<,'Iv(f^PFaPF«p)'Iv(f^PF^PF,p) 

+ ^ c^-CApC«<,c^^ TV (f^PFaPF^PF^p) Tr (f^PF^PF^PF^p) 

IJ.i'XpKaT^riC 

xTr (^F^PFjPr^PFpPr^p) 
- ^ '^^ Yl c^i^cxpC^TCr^Crjc TV (f^PFaPF^PF^p) TV (f.PF^PF^p) 

p,i'XpKaT^ri<^ 

xTV (r^vr^vr^v^ 

+ 0(A6), (A.6) 

where "TV " implies a trace taken with respect to the spinor indices. 

The complete list of 2PI planar vacuum diagrams at orders 6, 7, 8, 9 are given in Figs. 
7-13. The number of diagrams at each order is 4, 9, 24, 81, respectively. 

A.l Schwinger-Dyson equations 

In this section we derive a closed set of equations, which allows us to calculate J^{t, u) from 
G(c, d). 

Let be the one-particle-irreducible (IPI) part of the radiative corrections to c^i/. 
Then we can write c^^ as a geometric series 

111111 1 , 

c=- + -B- + -B-B- + --- = , (A.7) 
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Figure 6: The 2PI planar vacuum diagrams up to the 5-th order. For simphcity the propagators 
for bosons and fermions are represented by single (dashed and solid, respectively) lines. The first 
three diagrams arc contributions from the orders 1, 3, 4, and the last two from the order 5. (There 
is no contribution from the 2nd order.) The symmetry factor is indicated below each diagram. 
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Figure 7 : The 2PI planar vacuum diagrams at the 6-th order. 



where we regard c^,/, r^,y = t^S^i^, B^i, as 4 x 4 matrices. Note, on the other hand, that 
Bf^i^ can be obtained from G(c, d) by 



dc, 



Combining (A. 7) and (A. 8), we obtain 



d 



Similarly we get 



(2^~')a/3 = Aap + -^^G{c,d) , 



1 d 



which may be written in terms of as 



LA. 
2^dd^, 



G{c,d) , 



(A.8) 



(A.9) 



(A.IO) 



(A.ll) 



-13- 




Figure 8: The 2PI planar vacuum diagrams at tlie 7-th order. 



where A = E^('^m)^- 

Since the SD equations (A. 9), (A. 10) are closed with respect to c^y, T>a,/3, we can 
solve them order by order in A once the sum of the 2PI diagrams G(c, d) is obtained. For 
example, the results for the first few orders are as follows. 

oo 

c^. = 5] A" c(r.) , (A.12) 

n=0 

c(°) - - (5 

4i^ = -^ffTr(r^^-ir,^-i) , 

''U ''1/ ^ ' 



c(2) 


1 1 1 ... 
= -r > Tr 

p (-/i 

> Tr 

iff 
p '■M "-P "-i^ 










oo 
n=0 






(A.13) 


p(0) 


' J a/3 


5 
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Let us note that the free energy J^{t, u) satisfies the differential equations 

d 



2r du 



■r{t,u) 



(A.14) 
(A.15) 



By integrating these equations term by term, we finahy obtain the A-expansion for J-{t^ u). 
The integration constant, which appears only in the zcroth order contributions (i.e., the 
one-loop diagrams), can be fixed easily by direct computation. 
The explicit form of the free energy up to the 2nd order is 



A:=0 

2(1 - r) InA^ - 21n2 + ^ ^Ini^ - rln(det^) , 



J^{t, u) 

^0 



(A.16) 



:f2 



1 1 



A. 2 Exploiting the Ansatze 

As explained in Section 3 we impose some Ansatz on the possible symmetry breaking 
pattern in order to reduce the number of free parameters in the gaussian action. In this 
section we write down the SD equations for each Ansatz separately. By solving the SD 
equations, we obtain the free energy similarly to what we have done in section A.l. 

A. 2.1 SD equations for the SO(3) Ansatz 

The full propagators take the following form. 

/cn \ 

c 



\ 



do = d-i = cLa = Q 



(A.17) 



The SD equations (A. 9), (A. 11) are rewritten as 



1 

Cll 

1 



1 

dl 



d 



ti-2^[G(c,d)]so(3) , 
^^[G(c,d)]so(3) > 



dcii 
2d_ 

1 d 

+ — — [G(c,d)]so(3) • 



2r ddi 



(A.18) 
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Here and henceforth the symbol [ • ]so(n) iniphes that the number of independent variables 
is already reduced by the SO(n) Ansatz. The differential equations (A. 14), (A. 15) become 

2— [j^(t,u)]so(3) 
2d_ 

1 d 



Cll 



^[^{t^ ^)]S0(3) 



2r dui 

A.2.2 SD equations for the SO (2) Ansatz 

The full propagators take the following form. 

/ Cll Ci2 \ 
Cl2 Cll 



[•?^(*,^^)]so(3) = di ■ 



(A.19) 



V 



di = d2, ds = = . 



(A.20) 



7 



The SD equations (A. 9), (A. 11) are rewritten as 



Cll 



(cii)2 - (Ci2)2 
Cl2 



d 



dc 



11 



[G(c,d)]so(2) , 



(cii)' 



(C12)2 

1 



dc 



12 



[G(c, d)]so(2) , 



- = t 



1 

di 



— [G(c,d)]so(2) , 
1 d 

+ 2;:^[<^(c,t^)]so(2) 



(A.21) 



The differential equations (A. 14), (A. 15) become 



Cll 



at 
d 



r[-^(*'^)]so(2) = C 



4r dui 



[jr(t,ix)]go(2) = di . 



(A.22) 
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Figure 10: The 2PI planar vacuum diagrams at the 9-th order (to be continued). 
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Figure 11: The 2PI planar vacuum diagrams at the 9-th order (to be continued). 
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Figure 12: The 2PI planar vacuum diagrams at the 9-th order (to be continued). 
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Figure 13: The 2PI planar vacuum diagrams at the 9-th order. 



-23- 



